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Introduction 


Given  the  number  of  women  who  require  a  second  breast  cancer  removal  operation  due  to 
tumor-positive  margins,  and  the  limitations  of  techniques  currently  available  to  examine 
margins,  there  is  a  strong  need  for  a  device  that  can  quickly,  accurately,  and  non-invasively 
evaluate  the  margins  of  the  resected  breast  specimen  while  the  patient  is  in  the  operating  room. 
In  this  proposal,  we  hypothesize  that  spatially  offset  Raman  spectroscopy  (SORS)  can  be  used 
to  assess  the  margin  status  of  a  resected  breast  specimen  in  real  time,  in  the  operating  room 
with  high  sensitivity  and  specificity.  The  main  objectives  of  this  proposal,  then,  are  to 
characterize  the  depth  sensitivity  of  SORS  in  soft  tissues  and  implement  (i.e.  evaluate,  design 
and  test  this  technique  for  breast  tissues.  To  achieve  these  objectives,  the  following  specific 
aims  are  proposed:  (1)  Characterize  SORS  in  soft  tissues  -  the  relationship  between  source- 
detector  separation  and  depth  of  i  nterrogation  will  be  characterized  using  both  phantoms  and 
breast  tissue  samples.  (2)  Evaluate  the  performance  of  SORS  in  normal  and  malignant  breast 
tissues  in  vitro  -  Using  optical  fibers  and  the  results  of  aim  1,  the  ability  of  SORS  to  identify 
normal  and  malignant  breast  tissues  will  be  evaluated  and  compared  to  that  of  conventional  RS. 
(3)  Design,  build  and  test  a  1st  generation  SORS  device  on  resected  breast  specimens  in  the 
operating  room  -  Using  the  specifications  evaluated  under  aims  1  &  2,  a  SORS  fiber  device  will 
be  constructed  and  30  br  east  specimens  resected  from  partial  mastectomy  patients  will  be 
studied  at  multiple  depths  and  compared  with  histology. 
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Annual  Report 
Year  2 

(January  2010  -  December  2010) 


Body  of  the  report 

An  automated  method  that  can  assess  the  entire  resected  specimen  during  breast  conserving 
therapy  for  clear  margins  within  1-2  mm  of  the  surface  with  high  sensitivity  and  specificity  will 
assure  complete  removal  in  a  single  procedure  and  significantly  improve  management  of  breast 
cancer  therapy.  In  thi  s  grant,  we  proposes  to  develop  spatially  offset  Raman  spectroscopy 
(SORS)  to  detect  tumor  signatures  up  to  1  -2  mm  deep  and  thus  provide  real  time  feedback  to 
the  surgeon  so  that  the  complete  tumor  is  removed  in  a  s  ingle  visit.  Hence  we  propose  to 
develop  SORS  for  use  in  evaluating  surgical  margin  status  of  partial  mastectomy  specimens  in 
real  time  in  the  operating  room. 

In  year  2,  the  approved  statement  of  work  sought  to  develop  a  SORS  fiber  probe  for  breast 
tissues  and  to  develop  an  algorithm  that  optimizes  detection  of  positive  tumor  margins. 

Aim  2  -  Evaluate  the  performance  of  a  dedicated  SORS  probe  in  normal  and  malignant 
breast  tissues  in  vitro.  (Months  11-24). 

Task  3  -  Build  and  test  a  SORS  device 

1 .  Design  a  SORS  fiber  probe  with  the  source-detector  fiber  separations  identified  in  Aim  1 
and  other  parameters  such  as  fiber  diameter,  number  of  fibers  etc.  as  defined  above. 
(Months  11-12) 

2.  Build  the  probe  using  a  1 -dimensional  array  approach.  (Months  12-13) 

3.  Integrate  the  probe  with  a  diode  laser  and  Raman  spectrometer.  Timing  electronics  will 
be  developed  to  controlled  selection  of  fiber  combinations.  (Months  13-14) 

4.  Validate  the  SORS  system  using  the  ti  ssue  model  as  well  as  in  vitro  breast  tissue 
samples  for  performance  based  on  the  expected  projections  in  Aim  1.  Signal  strength, 
integration  time  and  processing  will  all  be  analyzed  for  optimal  output.  (Month  15) 

Task  4  -  Discrimination  and  Performance 

5.  Acquire  15-20  normal  and  malignant  breast  tissue  samples.  In  particular,  obtain  sample 
that  have  both  normal  and  malignant  regions  within  the  same  sample.  (Months  11-15) 

6.  Use  various  offsets  from  an  interrogation  depth  of  1  -2  mm  and  collect  spectra  from 
uniform  normal  and  malignant  breast  tissues  with  overlying  fat  or  normal  fibroglandular 
layers.  (Months  16-17) 

7.  Process  and  analyze  spectra  to  evaluate  the  integration  times,  signal  to  noise  and  depth 
response  from  normal  and  malignant  samples.  (Months18-20) 
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8.  Develop  multivariate  statistical  methods  using  spectra  from  1-2  mm  depths  for 
discrimination  between  normal  and  malignant  tissues.  Performance  estimates  will  be 
obtained  and  compared  to  previous  results.  (Months  21-22) 

9.  Investigate  the  applicability  of  the  technique  in  intact  heterogeneous  breast  tissue 
samples  that  have  normal  and  malignant  regions  using  the  SORS  device.  (Months  22) 

10.  Test  the  discrimination  algorithms  developed  above  to  obtain  the  performance  estimates 
of  SORS  with  these  samples.  Algorithms  will  be  modified  as  needed.  (Months  22-23) 

11.  Evaluate  accuracy  and  spatial  profile  of  spectral  results  by  comparing  depth  s  pecific 
spectroscopic  results  with  detailed  histology  at  each  depth.  (Months  23-24) 

12.  The  SORS  device  will  be  refined  as  needed  to  i  mprove  signal  collection  and  reduce 
integration  time.  Additional  Monte  Carlo  simulations  will  be  run  as  needed  to  ensure  that 
the  1-D  design  concept  can  be  adapted  to  interrogate  an  entire  surface  of  a  specimen 
with  similar  signal  to  noise  ratios.  (Months  22-24) 

Aim  2  Task  3:  The  primary  criterion  for  designing  a  SORS  probe  for  breast  tumor  margin 
analysis  was  to  ensure  proper  depth  sampling  -  that  is,  to  develop  a  probe  sensitive  to  tumor 
spectral  signatures  if  the  tumor  is  anywhere  within  the  fi  rst  2mm  in  depth  from  the  excised 
surface.  As  noted,  the  relevant  S-D  offsets  for  this  purpose  were  determined  to  be  <  3.5  mm. 
To  investigate  the  drop  in  SNR  as  S-D  offset  increases,  SORS  Monte  Carlo  simulations  were 
run  using  the  same  model  that  was  described  in  last  year’s  report  [1],  The  Raman  MC  code  is  a 
modified  version  of  a  multi-layer,  multi-fluorophore  code  developed  by  Vishwanath  and  Mycek 
[2-4].  The  code  uses  a  typical  Monte  Carlo  implementation  to  simulate  photon  propagation  and 
detection.  It  also  defines  a  probability  that  a  given  excitation  photon  will  be  Raman  scattered, 
according  to  the  tissue  layer's  properties,  at  each  propagation  step.  Simulations  were  run  for  3- 
layered  samples,  consisting  of  a  top  layer  of  0.5,  1 ,  or  2  mm  of  normal  breast  tissue,  a  0.1  to  20 

mm  thick  middle  layer  of  breast  tumor,  and  then  a 
2  cm  thick  bottom  layer  of  normal  breast  tissue  to 
mimic  the  clinical  situation  of  semi-infinite 
geometry.  Figure  1  shows  the  mean  of  signal  to 
noise  curves  obtained  using  the  3 -layer  model. 
As  predicted,  the  number  of  R  aman  photons 
detected  fell  off  at  w  hat  appears  to  be  an 
exponential  decay  as  a  function  of  S-D  offset 
regardless  of  relative  thicknesses  of  normal  and 
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cancerous  breast  tissues. 


These  results  were  used  to  design  a  multi-separation  SORS 
probe  (assembled  by  EMVision,  Loxahatchee,  FL),  whose 
distal  tip  is  shown  in  Figure  2.  A  single  400  pm  diameter 
source,  or  excitation  fiber  is  found  on  one  end,  and  f  our 
(partial)  rings  of  300  pm  diameter  collection  fibers  extend 
radially  outward.  The  excitation  fiber  includes  a  bandpas  s 
filter  at  its  tip  to  narrow  the  laser  line,  and  the  collection  fibers 
have  longpass  filters  at  their  tips  to  reject  elastically  scattered 
light.  The  center  to  center  distances  of  the  excitation  fiber  to 
each  detection  ring  are  0.5,  1.5,  2.5,  and  3.5  mm.  Based  on  Figure  1,  an  additional 
collection  fiber  was  added  to  each  consecutive  detector  ring  to  make  the  SNRs  from  each 
ring  more  comparable  to  one  another.  Adding  a  fiber  in  each  successive  detector  ring  had 
an  added  benefit  of  increasing  the  lateral  sampling  volume  of  the  probe  as  well. 

The  resulting  SORS  probe  was  connected  to  a  785  nm  diode  laser  (I0785MB0350M, 
Innovative  Photonics  Solutions,  Monmouth  Junction,  NJ)  which  delivered  80  mW  of  power 
to  the  sample.  The  collection  fibers  delivered  light  to  a  NIR-optimized  spectrograph  (LS785, 
Princeton  Instruments,  Princeton,  NJ),  which  dispersed  the  light  to  be  recorded  by  a  deep 
depletion,  thermo-electrically  cooled  CCD  (Pixis  400BR,  Princeton  Instruments).  A  neon- 
argon  lamp,  naphthalene,  and  acetaminophen  standards  were  used  to  c  alibrate  the 
wavenumber  axis,  and  a  N  IST-calibrated  tungsten-halogen  lamp  was  used  to  correct  for 
the  system  response  [5],  To  c  reate  a 
composite  spectrum  with  equal  weighting 
from  all  four  rings,  which  would  contain 
information  from  the  entire  2  mm 
sampling  depth,  the  binned  spectra  from 
each  of  the  four  rings  were  averaged 
after  processing. 


To  ensure  the  probe's  ability  to  gather 
spectra  from  each  ring  with  comparable 
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Figure  3:  Experimental  validation  of  SNR  from  SORS 
probe  in  chicken  breast. 


8 


SNRs,  spectra  were  acquired  for  20  seconds  each  from  12  different  spots  on  a  ~1  cm  thick 
piece  of  chicken  breast  (muscle).  Figure  3  shows  the  results  of  the  SNR  testing  on  chicken 
muscle.  Rings  1  and  4  of  the  SORS  probe,  with  one  and  four  fibers  per  ring,  and  with  S-D 
offsets  of  0.5  and  3.5  mm,  respectively,  displayed  nearly  identical  SNRs.  Rings  2  and  3 
showed  smaller  SNRs  compared  with  Ring  1,  but  only  by  -30%  and  20%,  respectively. 
This  trend  was  expected  based  on  the  shape  of  Figure  1,  indicating  that  the  design  of  the 
SORS  probe  effectively  accounted  for  SNR  fall-off  with  increasing  S-D  offset. 

Next  the  probe  was  tested  on  pure  normal  and  malignant  breast  tissue  samples  obtained 
with  required  approvals.  Figure  4  shows  typical  composite  spectra  recorded  from  pure 
normal  breast  tissue  and  a  pure  breast  tumor  (invasive  ductal  carcinoma)  with  the  SORS 
probe.  Numerous  spectral  regions  with  major  differences  between  the  two  tissue  types 
were  observed  and  proved  to  be  similar  to  the  spectra  reported  in  the  previous  annual 
report.  In  particular,  tumor  tissue  contains  a  strong  band  at  1006  cm'1,  usually  attributed  to 
phenylalanine,  while  normal  tissue  does  not.  The  r  atios  of  the  1  303  cm'1  to  1265  cm'1 
bands,  indicative  of  the  ratio  of  lipid  to 
protein  content,  are  very  different 
between  the  tissue  types,  and  the  amide 
I  band  c  entered  around  1656  cm'1  is 
much  wider  in  tumor  compared  to  normal 
-  again  indicative  of  i  ncreased  relative 
protein  contributions  in  the  c  ancerous 
tissues.  Also,  the  1445  cm'1  CH  stretch 
band  is  relatively  more  intense  in  normal 
tissue,  and  the  normal  tissue  contains  a 
carbonyl  stretch  peak  around  1750  cm"1, 
typically  due  to  f  at  content,  while  the 
tumor  tissue  does  not. 

Aim  2  Task  4:  In  total,  35  in  vitro  samples  from  total  mastectomy  patients  were  included  in 
the  study;  15  samples  had  either  no  tumor  (n=13)  or  tumor  >  2  mm  beneath  normal  tissue 
at  the  point  of  measurement  (n=2)  and  were  thus  labeled  as  "negative  margins,"  while  20 
samples  had  tumor  regions  within  the  first  2  mm  from  the  measurement  surface,  and  were 
thus  labeled  "positive  margins."  Wherever  possible,  measurements  from  tumor  samples 


9 


were  taken  such  that  the  SORS  probe  was  placed  on  a  small  region  of  visually  normal¬ 
appearing  tissue  on  top  oftheac  tual  tumor  to  mimic  the  situation  of  margin  evaluation. 
Spectra  were  recorded  for  10-30  seconds  and  processed.  Measurement  sites  were  inked, 
fixed  in  formalin,  and  serially  sectioned  to  correlate  the  spectra  with  histopathology 
diagnoses  of  tissue  type  and  the  depths  from  the  measurement  surface  of  those  tissues.  In 
this  manner,  the  analysis  was  done  to  di  scriminate  “negative”  margins  from  “positive” 
margins. 


Figure  5  shows  H&E  stained  tissue  section  and  the  SORS  spectra  from  a  sample  in  vitro  tumor 
sample.  The  "S"  arrow  in  the  histology  section  indicates  the  placement  of  the  source  fiber,  while 
the  "R1,"  "R2,"  etc.  labels  denote  the  location  of  the  individual  collection  fiber  rings.  The  sample 
provides  an  i  llustration  of  what  happens  with  smaller  layers  of  normal  tissue  over  a  turn  or. 
Figure  5A  shows  a  sample  with  a  large  tumor  region  but  with  pockets  of  normal  adipose  cells 
near  the  surface,  including  directly  under  the  location  where  the  excitation  fiber  from  the  probe 

was  placed.  Fr  om  Figure  5B,  the 
spectrum  from  the  smallest  S-D  offset 
contains  mostly  features  indicative  of 
normal  fatty  breast  tissue,  while  spectra 
from  the  I  arger  S-D  offsets  contain 
features  indicative  of  turn  or  spectral 
signatures.  Thus,  it  is  clear  that  the 
different  detector  rings  are  sampling 
different  volumes,  as  desired. 

To  simplify  the  “margin  analysis”  procedure,  the  spectra  from  each  detector  ring  were  averaged 
to  create  one  composite  spectrum  per  in  vitro  sample.  Table  1  shows  the  confusion  matrix  for 
classification  of  these  composite  spectra  with  Sparse  Multinomial  Linear  Regression,  a  non¬ 
linear  method  of  classification  developed  in  our  lab.  This  analysis  showed  an  excellent  ability 

for  SORS  to  evaluate  margin  status  in  breast 
specimens,  with  95%  sensitivity  and  100% 
specificity,  and  an  area  under  the  ROC  curve 
of  0.993.  Alternatively,  the  discrimination  was 
performed  with  94%  negative  predictive  value 
(NPV)  and  100%  positive  predictive  value 


Spectral  Margin  Status 

Negative 

Positive 

Histopathology 
Margin  Status 

Negative 

15 

1 

Positive 

0 

19 

Specificity:  100% 
Sensitivity:  95% 


NPV:  94%  PPV:  100% 

Table  1:  Confusion  matrix  showing  ability  of  SORS  to 
discriminate  positive  and  negative  margins. 
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(PPV).  This  initial  work  demonstrated  the  feasibility  and  promise  of  using  SORS  to  evaluate 
margin  status  on  intact  breast  specimens  in  a  laboratory  setting. 

Aim  2  Task  5  (NEW):  While  not  originally  proposed  in  the  grant  application  or  in  the  SOW,  a 
study  was  designed  to  test  the  feasibility  of  SORS  in  a  clinical  setting  using  the  same  SORS 
designed.  The  original  SOW  proposed  to  directly  design  a  2-D  surface  scanning  probe  before 
testing  SORS  in  a  clinical  intra-operative  setting.  However,  given  the  challenges  of  clinical 
application,  it  was  determined  thattesting  the  successful  1-D  probe  design  was  important  in 
developing  a  design  for  a  SORS  probe  that  would  be  able  to  scan  the  surface  of  breast  tissue. 
Thus,  a  30  patient  study  was  designed  where  SORS  measurements  were  performed  on  the 
surface  of  freshly  excised  partial  mastectomy  specimens  in  the  operating  room  in  real  time  while 
the  patient  was  still  present.  The  protocol  was  designed  such  that  one  measurement  was  made 
on  each  surgical  margin  (lateral,  superior,  deep,  etc.)  and  the  measurement  site  inked  for 
pathology.  These  inked  spots  were  subsequently  analyzed  using  histopathology  to  determine 
margin  status  of  these  areas  so  that  the  ability  of  SORS  to  detect  margin  may  be  evaluated. 


Figure  6  shows  the  composite  spectra  from 
a  lumpectomy  specimen  from  each  of  the  6 
sites  measured  as  shown  in  the  i  nset 
schematic.  The  spectra  inform  this  sample 
indicate  that  sites  1 ,  3,  and  5  are  from 
negative  margins  (location  5  is  from  a  very 
fibrous  area  with  scar)  and  sites  2,  4,  and  6 
are  from  margins  positive  for  DCIS.  A  total  of 
8  lumpectomies  were  included  in  this  study 
by  the  end  of  year  2  of  the  project. 


It  should  be  noted  that  at  this  time  the  graduate  student  working  on  the  project  (Matthew  Keller) 
completed  his  doctoral  thesis  and  graduated  from  Vanderbilt  University.  This  resulted  in  a  delay 
in  project  milestones  until  a  new  graduate  student  and  post-doctoral  fellow  could  be  recruited 
and  trained  for  the  completion  of  this  study  and  the  remaining  tasks  in  the  project. 


c 

■o  4 
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Raman  shift  (cm  ) 


Figure  6:  SORS  spectra  from  a  lumpectomy  specimen 
measured  in  the  OR.  Inset  shows  the  measurement 
protocol. 
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Key  Research  Accomplishments 

•  A  fiber  based  SORS  probe  was  designed 

•  This  SORS  probe  was  successfully  built  and  tested  in  chicken  breast  and  human  breast 
tissue  samples 

•  Signal  to  noise  estimates  from  the  simulations  were  experimentally  validated  and 
therefore  performance  of  the  SORS  probe  verified 

•  Breast  tissue  spectra  were  obtained  and  differences  observed  were  consistent  with 
previous  results 

•  SORS  probe  was  used  to  evaluate  “margins”  in  in  vitro  breast  tissue  specimen 

•  A  discrimination  algorithm  based  on  SMLR  was  developed  and  95%  sensitivity  and 
100%  specificity  was  observed  in  a  small  group  of  sample. 

•  A  new  study  was  designed  to  test  the  SORS  probe  in  a  clinical  setting  and  a  protocol  for 
study  developed 

•  8  lumpectomy  samples  were  successfully  evaluated  in  the  OR 
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Reportable  outcomes 


Conference  presentations 

•  Mahadevan-Jansen  A,  Keller  MD,  Kelley  M,  Clinical  implementation  of  SORS  for 
breast  tumor  margin  assessment,  Annual  meeting  of  the  Federation  of  Analytical 
Chemistry  and  Spectroscopy  Society  (FACSS),  Raleigh  NC,  October  2010. 

•  Mahadevan-Jansen  A,  Screening,  Diagnosis  and  G  uidance:  Can  Raman 
spectroscopy  help?,  Cancer  Colloqium,  St.  Andrews  Charitable  Trust,  St  Andrews 
UK,  November  2010. 

•  Mahadevan-Jansen  A,  Keller  MD,  Master  DB,  Kelley  M,  Spatially  offset  Raman 
spectroscopy  for  subsurface  assessment  of  breast  margins,  The  International 
Chemical  Congress  of  Pacific  Basin  Societies  (Pacifichem  2010),  Honolulu  HI, 
December  2010. 

Peer  reviewed  publications 

•  Keller  MD,  Majumder  SK,  Kelley  MC,  Meszoely  IM,  Boulos  FI,  Olivares  GM, 
Mahadevan-Jansen  A,  Autofluorescence  and  diffuse  reflectance  spectroscopy 
and  spectral  imaging  for  breast  margin  analysis,  Lasers  in  Surgery  and  Medicine, 
42(1),  15-23,  2010. 
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Conclusions 


The  project  has  the  tasks  and  sub-tasks  proposed  for  Aim  2  as  described  in  the  Statement  of 
Work  for  year  2  of  the  project.  We  were  able  to  successfully  build  a  SORS  fiber  probe  and  test  it 
on  tissue  phantoms,  standards  as  well  as  breast  tissue  samples.  The  probe’s  ability  to  assess 
margin  status  was  tested  using  an  in  vitro  breast  tissue  study  and  SORS  signal  from  35 
samples  were  analyzed.  A  discrimination  algorithm  was  developed  and  tested  and  found  to  be 
95%  sensitive  in  detecting  a  positive  margin.  To  validate  the  ability  of  SORS  in  a  clinical  setting, 
a  protocol  was  designed  to  measure  from  lumpectomy  specimens  while  the  patient  was  still  in 
operating  room.  Eight  lumpectomies  were  studied  with  excellent  preliminary  results  indicating 
the  promise  of  SORS  in  a  clinical  environment. 


14 


References 


1.  M.  D.  Keller,  R.  H.  Wilson,  M.  A.  Mycek  and  A.  Mahadevan-Jansen,  "Monte  Carlo  model 
of  spatially  offset  Raman  spectroscopy  for  breast  tumor  margin  analysis,"  Appl 
Spectrosc  64(6),  607-614  (2010) 

2.  K.  Vishwanath  and  M.  A.  Mycek,  "Time-resolved  photon  migration  in  bi-layered  tissue 
models,"  Opt  Express  13(19),  7466-7482  (2005) 

3.  K.  Vishwanath  and  M.  A.  Mycek,  "Do  fluorescence  decays  remitted  from  tissues 
accurately  reflect  intrinsic  fluorophore  lifetimes?,"  Optics  letters  29(13),  1512-1514 
(2004) 

4.  K.  Vishwanath,  B.  Pogue  and  M.  A.  Mycek,  "Quantitative  fluorescence  lifetime 
spectroscopy  in  turbid  media:  comparison  of  theoretical,  experimental  and  computational 
methods,"  Physics  in  medicine  and  biology  47(18),  3387-3405  (2002) 

5.  C.  A.  Lieber  and  A  .  Mahadevan-Jansen,  "Automated  method  for  subtraction  of 
fluorescence  from  biological  Raman  spectra,"  Appl  Spectrosc  57(11),  1363-1367  (2003) 


15 


Appendices 


1.  Keller  MD,  Vargis  E,  Granja  ND,  Robert  TH,  Mycek  MA,  Kelley  MC,  Mahadevan-Jansen 
A,  Development  ofaS  patially  Offset  Raman  Spectroscopy  Probe  for  Breast  Tumor 
Surgical  Margin  Evaluation,  Journal  of  Biomedical  Optics  16(7):  077006-1-8,  DOI: 
10.1117/1.3600708  (2011). 

Note:  This  paper  was  not  published  till  July  201 1 


16 


SPIEDigitalLibrary.org/jbo 


Development  of  a  spatially  offset  Raman 
spectroscopy  probe  for  breast  tumor 
surgical  margin  evaluation 


Matthew  D.  Keller 
Elizabeth  Yargis 
Nara  de  Matos  Granja 
Robert  H.  Wilson 
Mary-Ann  Mycek 
Mark  C.  Kelley 
Anita  Mahadevan- Jansen 


SPIE 


Downloaded  from  SPIE  Digital  Library  on  14  Feb  2012  to  129.59.95.1 15.  Terms  of  Use:  http://spiedl.org/terms 


Journal  of  Biomedical  Optics  1 6(7),  077006  (July  201 1 ) 


Development  of  a  spatially  offset  Raman  spectroscopy 
probe  for  breast  tumor  surgical  margin  evaluation 

Matthew  D.  Keller,3  Elizabeth  Vargis,3  Nara  de  Matos  Granja,b  Robert  H.  Wilson,0  Mary-Ann  Mycek,°'d'e  Mark  C.  Kelley/ 
and  Anita  Mahadevan-Jansen3 

aVanderbilt  University,  Department  of  Biomedical  Engineering,  Nashville,  Tennessee  37235 
bVanderbilt  University  Medical  Center,  Division  of  Pathology,  Nashville,  Tennessee  37232 
cUniversity  of  Michigan,  Applied  Physics  Program,  Ann  Arbor,  Michigan  48109 
^University  of  Michigan,  Department  of  Biomedical  Engineering,  Ann  Arbor,  Michigan  48109 
eUniversity  of  Michigan,  Comprehensive  Cancer  Center,  Ann  Arbor,  Michigan  48109 
fVanderbilt  University  Medical  Center,  Division  of  Surgical  Oncology,  Nashville,  Tennessee  37232 


Abstract.  The  risk  of  local  recurrence  for  breast  cancers  is  strongly  correlated  with  the  presence  of  a  tumor 
within  1  to  2  mm  of  the  surgical  margin  on  the  excised  specimen.  Previous  experimental  and  theoretical  results 
suggest  that  spatially  offset  Raman  spectroscopy  (SORS)  holds  much  promise  for  intraoperative  margin  analysis. 
Based  on  simulation  predictions  for  signal-to-noise  ratio  differences  among  varying  spatial  offsets,  a  SORS  probe 
with  multiple  source-detector  offsets  was  designed  and  tested.  It  was  then  employed  to  acquire  spectra  from  35 
frozen-thawed  breast  tissue  samples  in  vitro.  Spectra  from  each  detector  ring  were  averaged  to  create  a  composite 
spectrum  with  biochemical  information  covering  the  entire  range  from  the  tissue  surface  to  ~2  mm  below  the 
surface,  and  a  probabilistic  classification  scheme  was  used  to  classify  these  composite  spectra  as  "negative"  or 
"positive"  margins.  This  discrimination  was  performed  with  95%  sensitivity  and  100%  specificity,  or  with  100% 
positive  predictive  value  and  94%  negative  predictive  value.  © 2011  Society  of  Photo-Optical  Instrumentation  Engineers  (SPIE). 
[DOI:  10.1117/1.3600708] 

Keywords:  spatially  offset  Raman  spectroscopy;  breast  cancer;  margin  evaluation;  Monte  Carlo  simulation. 

Paper  1 1 041  R  received  Jan.  28,  2011;  revised  manuscript  received  May  6,  2011;  accepted  for  publication  May  23,  2011;  published 
online  Jul.  1,  201 1 . 


1  Introduction 

For  many  of  the  approximately  180,000  women  diagnosed  with 
early  stage  invasive  breast  cancer  or  carcinoma  in  situ  each  year, 1 
a  viable  treatment  option  is  breast  conserving  therapy  (BCT). 
The  surgical  portion  of  BCT  involves  a  partial  mastectomy,  or 
lumpectomy,  to  remove  only  the  primary  lesion  with  a  small 
amount  of  surrounding  normal  tissue.2  Depending  on  the  hospi¬ 
tal,  the  depth  of  normal  tissue  required  from  the  surgical  margin 
(i.e.,  the  surface  of  the  excised  specimen)  to  the  tumor  is  typi¬ 
cally  1  to  2  millimeters/  This  situation  is  illustrated  in  Fig.  1 ;  if 
a  sufficient  amount  of  normal  tissue  exists,  as  in  the  right  side  of 
Fig.  1,  the  margins  are  said  to  be  negative  for  tumor.  If  tumor¬ 
positive  margins  are  found,  as  shown  on  the  left  side  of  Fig.  1, 
a  second  operation  is  necessary  because  positive  margins  are  a 
major  predictor  of  local  tumor  recurrence.4  Currently  available 
intraoperative  margin  evaluation  tools,  such  as  simple  visual  ex¬ 
amination,  ultrasound,  cytological  examination  (“touch  prep”), 
and  frozen  section  analysis,  all  have  significant  drawbacks  in 
terms  of  accuracy  and/or  time  required,3-3  so  there  is  a  need  for 
an  automated,  real-time  method  to  accurately  evaluate  surgical 
margins  during  BCT. 

Several  optical  techniques  have  recently  been  investigated 
for  breast  tumor  margin  evaluation.  A  diffuse  reflectance  imag¬ 
ing  system9-11  was  employed  to  assess  55  tumor  margins  from 
48  patients,  and  classified  whole  margins  as  negative  versus 
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positive/close  with  79%  sensitivity  and  67%  specificity.11 
Kennedy  et  al.  employed  the  same  diffuse  reflectance  spec¬ 
troscopy  system  to  characterize  tumor  margins  from  100 
patients.12  Keller  et  al.  used  combined  diffuse  reflectance  and 
autofluorescence  spectroscopy  to  classify  individual  points  from 
32  patients  on  margins  as  negative  versus  positive  with  85%  sen¬ 
sitivity  and  96%  specificity,  and  also  demonstrated  the  ability  to 
perform  autofluorescence  spectral  imaging  of  larger  regions.13 
Nguyen  et  al.  used  optical  coherence  tomography  (OCT)  im¬ 
ages  from  20  patients  to  classify  the  margin  status  with  100% 
sensitivity  and  82%  specificity.14 

In  the  one  report  of  using  Raman  spectroscopy  for  a  margin 
analysis  tool,13  measurements  were  made  in  vivo  rather  than  on 
the  excised  specimen;  the  latter  approach  is  the  current  standard 
practice  in  surgical  pathology.  That  study  also  used  a  conven¬ 
tional  fiber  optic  probe  with  source  and  detector  fibers  adjacent 
to  each  other,  which  allows  only  limited  depth  sensing.  More 
recently,  spatially  offset  Raman  spectroscopy  (SORS)  has  been 
shown  to  be  a  reliable  method  for  recovering  biological  Raman 
spectra  from  depths  greater  than  those  possible  with  standard 
techniques. 16-24  It  does  so  because  detection  elements  spaced  ra¬ 
dially  further  from  source  elements  are  more  sensitive  to  photons 
traveling  deeper  beneath  the  tissue  surface  and  to  greater  radial 
distances  due  to  multiple  scattering  (demonstrated  in  Fig.  1). 
We  have  previously  demonstrated  the  ability  of  SORS  to  detect 
spectral  contributions  from  breast  tumors  buried  under  0.5  to 
2  mm  of  normal  breast  tissue.16  In  addition,  a  SORS  Monte 
Carlo  (MC)  code  was  developed  to  quantify  signals  obtained 
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Detectors  Detectors 


Fig.  1  Drawing  of  tumor-positive  versus  tumor-negative  margins,  de¬ 
fined  by  the  distance  between  the  surgical  margin  and  the  tumor 
boundary.  Overlaid  are  general  photon  migration  paths  demonstrat¬ 
ing  the  advantage  of  SORS  for  this  application. 


08 

Q  E 

8  = 

O  C 

o  x  0.6  - 

CL  5 
c  o 


0.4 


0.2  - 


1  2  3  4  5 

Source -Detector  Offset  (mm) 


Fig.  2  Simulation  results  for  total  number  of  Raman  photons  detected 
as  a  function  of  detector  S-D  offset,  normalized  to  a  maximum  of  1. 
Results  were  averaged  for  a  variety  of  thicknesses  of  the  top  two  tissue 
layers  (0.5  to  2  mm  for  top  normal  layer,  0.1  to  20  mm  for  middle 
tumor  layer).  No  error  bars  are  shown  because  the  standard  deviation 
was  less  than  1  %. 


from  layered  constructs  of  normal  breast  tissue  overlying  breast 
tumors.25  In  particular,  the  code  was  used  to  examine  the  effects 
of  layer  thicknesses  and  overall  geometries  on  relative  tumor 
contributions  to  detected  spectra  for  a  range  of  source-detector 
(S-D)  offsets.25  To  detect  a  tumor  signature  within  the  first  2 
mm  from  the  surface,  the  resulting  spectrum  at  a  given  S-D  off¬ 
set  must  contain  at  least  a  5%  contribution  from  the  tumor.  To 
achieve  this  level  of  contribution,  it  was  found  that  the  tumor 
would  have  to  be  ~0.1-mm  thick  under  0.5  mm  of  normal  tissue, 
or  ~l-mm  thick  under  2-mm  of  normal  tissue. 2:1 

Combining  the  results  of  experimental16  and  numerically 
simulated2"1  SORS  indicated  that  to  be  sensitive  to  breast  tumors 
located  up  to  2  mm  beneath  normal  breast  tissue,  as  needed  for 
surgical  margin  evaluation,  a  maximum  S-D  offset  of  ~3.5  mm 
should  be  used.  With  larger  offsets,  the  measurements  could 
possibly  detect  large  tumors  from  over  2  mm  deep  and  create 
false  positives;  also,  recording  spectra  with  adequate  signal- 
to-noise  ratios  (SNRs)  at  larger  offsets  is  difficult  since  fewer 
photons  tend  to  escape  the  tissue  surface  as  the  S-D  offset  in¬ 
creases.  In  addition,  a  shortcoming  of  the  previous  experimen¬ 
tal  work  was  the  need  to  translate  the  single  detector  fiber  for 
each  measurement.  Thus,  the  goal  of  this  work  was  to  design, 
test,  and  implement  a  multiseparation  SORS  probe  for  breast 
tumor  surgical  margin  evaluation.  In  particular,  this  manuscript 
describes  the  use  of  the  previously  developed  SORS  Monte 
Carlo  code  to  investigate  the  theoretical  drop  in  SNR  as  a  func¬ 
tion  of  the  S-D  offset,  the  design  of  a  SORS  probe  based  on  the 
above  theoretical  and  experimental  findings,  its  testing  to  ensure 
comparable  signal  quality  in  each  detector  ring,  and  its  use  in 
acquiring  spectra  from  breast  cancer  samples  to  assess  its  ability 
to  accurately  evaluate  surgical  margin  status. 

2  Materials  and  Methods 

2.1  SORS  Probe  Design 

The  primary  criterion  for  designing  a  SORS  probe  for  breast  tu¬ 
mor  margin  analysis  was  to  ensure  proper  depth  sampling — that 
is,  to  develop  a  probe  sensitive  to  tumor  spectral  signatures  if 


the  tumor  is  anywhere  within  the  first  2  mm  in  depth  from  the 
excised  surface.  As  noted,  the  relevant  S-D  offsets  for  this  pur¬ 
pose  were  determined  to  be  <3.5  mm.  To  investigate  the  drop 
in  SNR  as  S-D  offset  increases,  SORS  Monte  Carlo  simulations 
were  run  using  the  same  model  as  for  the  previous  results.2"1 
The  Raman  MC  code  is  a  modified  version  of  a  multilayer,  mul- 
tifluorophore  code  developed  by  Vishwanath  and  Mycek.26-28 
As  previously  described,2"1  the  code  uses  a  typical  Monte  Carlo 
implementation  to  simulate  photon  propagation  and  detection. 
It  also  defines  a  probability  that  a  given  excitation  photon  will 
be  Raman  scattered,  according  to  the  tissue  layer’s  properties,  at 
each  propagation  step.  Raman  scattered  photons  are  then  prop¬ 
agated  according  to  new  optical  properties  of  each  tissue  layer, 
and  those  Raman  photons  exiting  the  top  surface  of  the  tissue 
construct  within  the  collection  cone  of  a  given  detector  bin  are 
recorded.  Simulations  were  run  for  3-layered  samples,  consist¬ 
ing  of  a  top  layer  of  0.5,  1 ,  or  2  mm  of  normal  breast  tissue,  a  0. 1 
to  20-mm  thick  middle  layer  of  breast  tumor,  and  then  a  2-cm 
thick  bottom  layer  of  normal  breast  tissue  to  mimic  the  clinical 
situation  of  semi-infinite  geometry. 

As  a  metric  for  SNR,  the  total  number  of  simulated  Raman 
photons,  originating  from  any  layer,  reaching  each  detector  bin 
was  counted  and  normalized  to  a  maximum  of  1,  since  we  are 
only  interested  in  how  SNR  falls  off  with  the  S-D  offset.  Because 
the  raw  signal  strength  was  consistent  among  the  4  rings,  it 
was  assumed  that  noise  levels  were  consistent  among  the  4 
rings  as  well.  Figure  2  shows  the  mean  of  these  SNR  curves; 
since  the  standard  deviation  was  less  than  1%  over  the  range 
of  thicknesses  for  the  top  two  layers,  no  error  bars  are  shown. 
As  predicted,  the  number  of  Raman  photons  detected  fell  off  at 
what  appears  to  be  an  exponential  decay  as  a  function  of  the  S-D 
offset.  Although  the  trend  of  SNR  as  a  function  of  S-D  offset 
was  consistent  regardless  of  relative  thicknesses  of  normal  and 
cancerous  breast  tissues,  the  trend  may  not  hold  for  other  tissue 
types,  especially  those  in  which  optical  properties  can  vary  more 
drastically  in  inhomogeneous  regions. 

These  results  were  used  to  aid  in  the  design  of  a  multisep¬ 
aration  SORS  probe  (assembled  by  EMVision,  Loxahatchee, 
Florida),  whose  distal  tip  is  shown  in  Fig.  3.  A  single  400-/zm 
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Fig.  3  Schematic  of  tip  of  SORS  probe.  S — source  fiber,  all  other 
circles — collection  fibers.  Each  ring  of  collection  fibers  is  separated  by 
1  mm  (center  to  center)  from  the  previous  ring,  with  0.5  mm  between 
the  centers  of  the  source  fiber  and  the  first  detector  ring  (R1 ). 

diameter  source,  or  excitation  fiber,  is  found  on  one  end  and  4 
(partial)  rings  of  300-/xm  diameter  collection  fibers  extend  radi¬ 
ally  outward.  The  excitation  fiber  includes  a  bandpass  filter  at  its 
tip  to  narrow  the  laser  line,  and  the  collection  fibers  have  long- 
pass  filters  at  their  tips  to  reject  elastically  scattered  light.  The 
center  to  center  distances  of  the  excitation  fiber  to  each  detection 
ring  are  0.5,  1.5,  2.5,  and  3.5  mm.  Based  on  Fig.  2,  an  additional 
collection  fiber  was  added  to  each  consecutive  detector  ring  to 
make  the  SNRs  from  each  ring  more  comparable  to  one  another. 
While  the  curve  in  Fig.  2  is  not  linear,  adding  a  single  fiber 
for  each  larger-offset  ring  provided  the  closest  approximation  of 
equivalent  SNRs,  if  all  of  the  fibers  in  each  ring  were  binned. 
Adding  a  fiber  in  each  successive  detector  ring  had  an  added 
benefit  of  increasing  the  lateral  sampling  volume  of  the  probe 
as  well. 

2.2  Instrumentation  and  Data  Processing 

The  SORS  probe  delivered  80  mW  of  power  from  a  785  nm 
diode  laser  (I0785MB0350M,  Innovative  Photonics  Solutions, 
Monmouth  Junction,  New  Jersey).  The  collection  fibers  deliv¬ 
ered  light  to  a  near-infrared-optimized  spectrograph  (LS785, 
Princeton  Instruments,  Princeton,  New  Jersey),  which  dispersed 
the  light  to  be  recorded  by  a  deep  depletion,  thermo-electrically 
cooled  CCD  (Pixis  400BR,  Princeton  Instruments). 

Each  acquisition  with  the  SORS  probe  recorded  4  spectra — 
one  from  each  detector  ring.  Each  ring  was  calibrated  separately 
since  the  inherent  curvature  in  the  detection  system  created 
slight  but  noticeable  differences  in  peak  locations  on  the  CCD 
among  different  rings.  A  neon-argon  lamp,  naphthalene,  and 
acetaminophen  standards  were  used  to  calibrate  the  wavenum¬ 
ber  axis,  and  a  National  Instrument  Standards  and  Tecnology- 
calibrated  tungsten-halogen  lamp  was  used  to  correct  the  sys¬ 
tem  response.24  After  wavenumber  binning  (in  3.5  cm-1  steps 
given  the  system  resolution  of  ~7  cm-1)  and  noise  smoothing, 
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the  background  fluorescence  was  subtracted  with  a  modified 
polynomial  fit  algorithm,24  and  the  spectra  were  normalized  ac¬ 
cording  to  their  overall  mean  intensities.  To  create  a  composite 
spectrum  with  equal  weighting  from  all  4  rings,  which  would 
contain  information  from  the  entire  2  mm  sampling  depth,  the 
binned  spectra  from  each  of  the  4  rings  were  averaged  after 
processing. 


2.3  SNR  Testing 

To  ensure  the  probe’s  ability  to  gather  spectra  from  each  ring 
with  comparable  SNRs,  spectra  were  acquired  for  20  s  each 
from  12  different  spots  on  a  ~1  cm  thick  piece  of  chicken 
breast  (muscle).  The  spectra  were  processed  as  described  in 
Sec.  2.2,  and  the  SNR  of  the  binned  spectrum  from  each  ring 
was  calculated  by  dividing  the  height  of  the  1445  cm-1  peak, 
which  is  the  strongest  peak  in  all  samples  measured,  by  the 
standard  deviation  of  the  flat  (i.e.,  no  Raman  signal)  spectral 
range  between  the  1656  and  1750  cm-1  peaks,  which  represent 
the  noise  inherent  in  the  system  that  could  not  be  removed  via 
pre-  or  post-processing. 

2.4  In  vitro  Sample  Measurements 

With  approval  by  the  Vanderbilt  Institutional  Review  Board  (No. 
050551)  and  the  U.S.  Army  Medical  Research  and  Materiel 
Command’s  Human  Research  Protection  Office,  fresh  frozen 
human  breast  tissue  samples  were  acquired  from  the  Cooper¬ 
ative  Human  Tissue  Network.  Frozen-thawed  tissues  are  not 
perfect  surrogates  for  freshly  excised  tissue,  as  their  optical 
properties  can  differ  from  each  other.30  A  recent  study  by  Reble 
et  al.31  demonstrated  that  Raman  sampling  volumes  can  vary 
substantially  based  on  a  tissue’s  optical  properties,  especially 
the  reduced  scattering  coefficient.  Nevertheless,  using  such  tis¬ 
sues  is  a  common  first  step  for  breast  cancer  studies,30,32  and  a 
recent  study  showed  nearly  an  equivalent  performance  of  an  al¬ 
gorithm  for  differentiating  normal,  benign,  and  malignant  breast 
tissues  developed  with  in  vitro  Raman  spectra  and  applied  to  in 
vivo  measurements.33 

In  total,  35  samples  were  included  in  the  study;  15  samples 
had  either  no  tumor  ( n  =  13)  or  a  tumor  >  2  mm  beneath  normal 
tissue  at  the  point  of  measurement  (n  =  2),  and  were  thus  labeled 
as  “negative  margins,”  while  20  samples  had  a  tumor  [15  inva¬ 
sive  ductal  carcinoma  (IDC)  and  5  invasive  lobular  carcinoma 
(ILC)]  regions  within  the  first  2  mm  from  the  measurement 
surface,  and  were  thus  labeled  “positive  margins.”  Of  the  nega¬ 
tive  samples,  7  were  predominantly  adipose  and  8  were  varying 
compositions  of  adipose  and  fibroglandular  tissue.  Of  the  posi¬ 
tive  samples,  8  had  tumors  underlying  various  compositions  of 
normal  tissue  ranging  from  0.1  to  1.5  mm  thick,  and  12  samples 
had  tumor  regions  at  the  surface  under  at  least  part  of  the  probe. 
Wherever  possible,  measurements  from  the  tumor  samples  were 
taken  such  that  the  SORS  probe  was  placed  on  a  small  region  of 
visually  normal-appearing  tissue  on  top  of  the  actual  tumor  to 
mimic  the  situation  of  margin  evaluation.  Spectra  were  recorded 
for  10  to  30  s  and  processed  as  described  above.  Measurement 
sites  were  inked,  fixed  in  formalin,  and  serially  sectioned  to  cor¬ 
relate  the  spectra  with  histopathology  diagnoses  of  tissue  type 
and  the  depths  from  the  measurement  surface  of  those  tissues. 
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In  this  manner,  the  analysis  was  done  to  discriminate  “negative” 
margins  from  “positive”  margins. 

2.5  Classification  of  Margin  Status 

The  composite  spectrum  from  averaging  all  4  detector  rings 
was  used  for  analysis,  and  if  there  were  histological  evidence  of 
tumor  cells  within  2  mm  of  the  measurement  surface,  the  “mar¬ 
gin”  was  considered  positive.  All  tumor-positive  measurements 
were  lumped  into  a  single  category  since  the  surgeon  simply 
needs  to  know  whether  any  types  of  malignant  cells  remain  too 
close  to  the  margin.  Discrimination  was  performed  with  sparse 
multinomial  logistic  regression  (SMLR),34  a  Bayesian  machine¬ 
learning  framework  that  computes  the  posterior  probability  of 
a  spectrum  belonging  to  each  tissue  class  based  on  a  labeled 
training  set.  In  the  case  of  this  binary  analysis,  whichever  class 
had  the  higher  probability  of  membership  was  the  one  to  which 
the  spectrum  was  classified.  SMLR  also  includes  inherent  di¬ 
mensionality  reduction  as  it  seeks  to  create  sparse  basis  vectors, 
which  is  important  for  these  data  sets  given  their  small  sizes. 
Since  each  in  vitro  sample  had  only  one  measurement  site  (their 
sizes  relative  to  the  probe  precluded  multiple  independent  sites), 
SMLR  was  run  with  leave-one-out  cross-validation.  A  range  of 
input  parameters  to  SMLR  were  tested,  and  the  combination  that 
provided  the  most  accurate  classification,  while  also  maximiz¬ 
ing  sparsity,  was  using  a  Laplacian  prior,  direct  kernel,  lambda 
value  of  0.01,  and  not  adding  a  bias  term. 

3  Results 

Figure  4  shows  the  results  of  the  SNR  testing  on  the  chicken 
muscle.  Rings  1  and  4  of  the  SORS  probe,  with  1  and  4  fibers 
per  ring,  and  with  S-D  offsets  of  0.5  and  3.5  mm,  respectively, 
displayed  nearly  identical  SNRs.  Rings  2  and  3  showed  smaller 
SNRs  compared  with  ring  1,  but  only  by  ~30%  and  20%,  re¬ 
spectively.  This  trend  was  expected  based  on  the  shape  of  Fig.  2, 
although  the  signal  strengths  of  rings  2  and  3  were  smaller  than 
predicted  by  the  simulations.  The  likely  reason  is  that  when 
imaging  the  detection  fibers  for  the  two  middle  rings  during 
alignment  and  testing,  their  throughput  appeared  to  be  lower, 
compared  with  the  fibers  for  rings  1  and  4.  Even  so,  the  design 
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Fig.  4  Mean  (n  =  12)  SNR,  plus  or  minus  one  standard  deviation, 
for  spectra  of  chicken  muscle  binned  within  each  detector  ring  and 
normalized  to  SNR  of  first  ring. 
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Fig.  5  Typical  normalized  composite  spectra  from  SORS  probe  of  nor¬ 
mal  breast  tissue  versus  malignant  breast  tumor  tissue. 


of  the  SORS  probe  effectively  accounted  for  SNR  fall-off  with 
increasing  S-D  offset. 

Figure  5  shows  typical  composite  spectra  recorded  from  pure 
normal  breast  tissue  and  a  pure  breast  tumor  (invasive  ductal  car¬ 
cinoma  for  this  example;  other  tumors  have  similar  spectra35) 
tissue  with  the  SORS  probe.  As  in  the  previous  study,1"  there 
are  numerous  spectral  regions  with  major  differences  between 
the  two  tissue  types.  In  particular,  tumor  tissue  contains  a  strong 
band  at  1006  cm-1,  usually  attributed  to  phenylalanine,  while 
normal  tissue  does  not.  The  ratios  of  the  1303  to  1265  cm-1 
bands,  indicative  of  the  ratio  of  lipid  to  protein  content,  are  very 
different  between  the  tissue  types,  and  the  amide  I  band  centered 
around  1656  cm-1  is  much  wider  in  the  tumor  as  compared  to 
normal — again  indicative  of  increased  relative  protein  contribu¬ 
tions  in  the  cancerous  tissues.  Also,  the  1445  cm-1  CH  stretch 
band  is  relatively  more  intense  in  normal  tissue,  and  the  nor¬ 
mal  tissue  contains  a  carbonyl  stretch  peak  around  1750  cm-1, 
typically  due  to  fat  content,  while  the  tumor  tissue  does  not. 

Figures  6  and  7  show  H&E  stained  tissue  sections  and  the 
SORS  spectra  from  those  sections  from  three  different  [Figs.  6, 
7(a)  and  7(b),  and  7(c)  and  7(d)]  in  vitro  tumor  samples.  In  all 
histological  images,  the  “.S’”  arrow  indicates  the  placement  of 
the  source  fiber,  while  the  “Rl,”  “R2,”  etc.,  labels  denote  the 
location  of  the  individual  collection  fiber  rings.  In  the  tissue 
sample  from  Fig.  6(a),  the  probe  was  delivering  light  to  a  large 
fatty  area,  as  seen  by  the  whitish  (formerly  lipid-filled)  vacuoles, 
while  only  the  outermost  collection  fibers  were  placed  over  a 
portion  of  the  tumor,  which  comprises  the  remainder  of  the 
darkly  stained  section.  Since  spectral  differences  among  detec¬ 
tor  rings  in  Fig.  6(b)  are  visually  subtle,  except  for  differences 
around  1445  cm-1,  close-ups  of  the  three  spectral  regions  are 
shown  in  Figs.  6(c)-6(e).  These  plots  show  definite  trends  in¬ 
dicating  that  the  closer  rings  are  sampling  normal  tissue,  while 
the  outer  rings  are  picking  up  slight  spectral  contributions  from 
the  tumor  as  well.  By  qualitatively  comparing  the  spectra  in 
Fig.  6(b)  with  the  typical  pure  normal  and  tumor  spectra  from 
Fig.  5,  these  trends  include  the  increasing  presence  of  the  1006 
cm-1  peak,  the  lesser  relative  contributions  from  the  1303  and 
1445  cm-1  peaks,  and  the  increasing  width  of  the  1656  cm-1 
peak  as  the  source-detector  offset  increases.  These  trends  are 

July  2011  •  Vol.  16(7) 


077006-4 


Downloaded  from  SPIE  Digital  Library  on  14  Feb  2012  to  129.59.95.1 15.  Terms  of  Use:  http://spiedl.org/terms 


Keller  et  al.:  Development  of  a  spatially  offset  Raman  spectroscopy  probe... 


(a) 


R4  R3  R2  R1  S 

—trill 


(b) 


Raman  shill  (cm'1) 


(C) 


(d) 


(e) 


Fig.  6  (a)  H&E  stained  tissue  section  of  IDC  sample  with  a  large  area  of  normal  fat  (white  colored  area  with  "N")  on  the  right,  and  solid  IDC  tumor 
(darkly  stained  area  with  "T")  on  the  left.  Arrows  indicate  the  placement  of  the  source  fiber  (S)  and  each  of  the  detector  rings,  (b)  Binned  SORS 
spectra  for  each  detector  ring  from  tissue  in  A.  To  aid  in  the  visualization  of  relevant,  but  subtle  spectral  changes,  zoomed-in  versions  are  shown 
for  (c)  the  1 006  cm  _1  phenylalanine  peak,  (d)  the  1265  cm-1  amide  III  and  1 303  cm-1  lipid  peak,  and  (e)  the  shoulders  of  the  1  656  cm amide 
I  peak. 


similar  to  those  seen  in  the  earlier  report  of  SORS  on  layered 
breast  tissues,16  but  in  this  case,  the  tissue  boundary  was  vertical 
rather  than  horizontal. 

The  example  in  Figs.  7(a)  and  7(b)  provides  an  illustration 
of  what  happens  with  smaller  layers  of  normal  tissue  over  a 
tumor.  Figure  7(a)  shows  a  sample  with  a  large  tumor  region,  but 
with  pockets  of  normal  adipose  cells  near  the  surface,  including 
directly  under  the  location  where  the  excitation  fiber  from  the 
probe  was  placed.  From  Fig.  7(b),  in  comparison  to  Fig.  5,  the 
spectrum  from  the  smallest  S-D  offset  mostly  contains  features 
indicative  of  normal  fatty  breast  tissue,  while  spectra  from  the 
larger  S-D  offsets  contain  features  indicative  of  tumor  spectral 
signatures,  as  noted  above.  The  sample  from  Figs.  7(c)  and 
7(d)  is  included  to  confirm  that  if  the  excitation  side  of  the 
probe  is  placed  on  the  tumor  tissue  overlying  normal  tissue  (i.e., 
the  opposite  of  margin  analysis),  then  the  inner  detector  rings 
picked  up  tumor  signatures,  while  the  outer  rings  picked  up  the 
appropriate  degree  of  normal  spectral  signatures.  Thus,  it  is  clear 
that  the  different  detector  rings  are  sampling  different  volumes, 
as  desired. 

To  simplify  the  “margin  analysis”  procedure,  the  spectra  from 
each  detector  ring  were  averaged  to  create  one  composite  spec¬ 
trum  per  in  vitro  sample.  Thus,  a  single  histological  classifi¬ 
cation  could  be  correlated  to  a  single  spectral  classification. 
Table  1  shows  the  confusion  matrix  for  the  classification  of 
these  composite  spectra  with  SMLR.  This  analysis  showed 
an  excellent  ability  for  SORS  to  evaluate  margin  status  in 
breast  specimens,  with  95%  sensitivity  and  100%  specificity, 
and  an  area  under  the  receiver  operating  characteristic  curve 
of  0.993.  Alternatively,  the  discrimination  was  performed  with 

Journal  of  Biomedical  Optics 


94%  negative  predictive  value  (NPV)  and  100%  positive  pre¬ 
dictive  value  (PPV).  The  one  false  negative  came  from  a  tu¬ 
mor  sample  which,  after  formalin  fixation  and  sectioning,  was 
found  to  have  a  ~1.5  mm  layer  of  normal  tissue  between 
the  measurement  site  and  the  tumor.  Since  it  has  been  shown 
that  normal  tissue  margins  tend  to  shrink  by  an  average  of 
33%  during  formalin  fixation,36  it  is  possible  that  this  nor¬ 
mal  layer  was  at  least  2-mm  thick  when  the  spectra  were 
obtained. 

4  Discussion 

This  manuscript  presents  the  design,  testing,  and  implemen¬ 
tation  of  a  multiseparation  SORS  probe  for  use  in  evaluating 
surgical  margin  status  following  partial  mastectomies.  The  de¬ 
sign  shown  in  Fig.  3  was  based  on  the  results  from  our  earlier 
experimental  and  simulation-based  studies,16-23  and  from  the 
SNR  simulation  results  from  Fig.  2.  To  ensure  that  the  SNRs 
were  comparable  across  the  different  detector  rings,  a  series  of 
measurements  was  performed  using  the  common  soft  tissue  op¬ 
tical  phantom  of  a  chicken  breast.  As  seen  in  Fig.  4,  the  design 
of  adding  an  additional  collection  fiber  for  each  further-offset 
ring  worked  well  to  keep  the  SNR  of  each  ring  no  more  than 
~30%  different  from  the  others.  Given  the  exponential  shape  of 
Fig.  2,  it  would  be  very  difficult  to  design  a  probe  to  both  sample 
the  desired  depths  in  tissue  and  achieve  even  better  equilibra¬ 
tion  of  SNR  among  the  various  detector  rings.  Besides  the  SNR 
balancing,  the  probe  design  from  Fig.  3  also  appeared  to  sample 
tissue  to  the  expected  depths  based  on  earlier  experimental16 
and  simulation25  results.  This  conclusion  is  supported  by  the 
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(a)  (b) 
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Fig.  7  (a)  H&E  stained  tissue  section  of  ILC  sample  with  pockets  of  normal  fat  ("N")  near  surface  of  otherwise  darkly  stained  tumor  tissue  ("T"). 
Arrows  again  represent  placements  of  fibers,  (b)  SORS  spectra  for  each  detector  ring  from  tissue  in  A.  (c)  H&E  stained  tissue  section  of  IDC  sample 
with  underlying  fat  and  (d)  corresponding  SORS  spectra. 


success  shown  in  Table  1  for  classifying  spectra  according  to 
the  margin  status  using  2  mm  as  the  cutoff  value  for  negative 
versus  positive  classification. 

The  ability  of  the  detector  rings  to  sample  different  volumes 
is  demonstrated  in  Figs.  6  and  7.  From  Fig.  6(a),  it  is  clear  that  the 
SORS  probe  was  placed  over  two  very  different  regions  of  tissue 
for  that  specimen.  A  large  area  of  normal  fatty  tissue  was  found 
directly  under  the  excitation  fiber  and  the  first  2  to  3  detector 
rings,  while  the  outermost  1  to  2  detector  rings  were  placed 
against  the  tumor.  Comparing  Figs.  6(b)-6(e)  to  the  pure  normal 
and  tumor  spectra  from  Fig.  5,  rings  1  and  2  show  essentially  no 
tumor  spectral  signatures.  Given  this,  a  standard  Raman  probe 
placed  in  the  same  spot  would  not  detect  any  positive  margin 


Table  1  Confusion  matrix  for  "margin  analysis"  on  in  vitro  specimens. 


Spectral  margin  status 

Negative  Positive 

Histopathology  Negative  15  0 

Specificity:  100% 

Margin  Status  Positive  1  19 

NPV:  94%  PPV:  100% 

Sensitivity:  95% 

findings  at  this  point.  The  3rd  and  4th  rings  of  the  probe  were 
able  to  detect  slight  tumor  contributions  though,  indicating  that 
they  successfully  sampled  a  different  volume  of  tissue  than  the 
inner  rings.  While  most  spectral  regions  showed  an  increasing 
tumor  contribution  from  ring  1  to  ring  4,  ring  3  had  a  stronger 
relative  1006  cm-1  peak  than  ring  4.  Possible  causes  for  this 
include  a  slight  misalignment  between  the  probe  and  histology, 
especially  considering  the  rotation  of  the  probe,  or  inconsistent 
biochemical  composition  of  the  tumor  tissue  sampled  by  rings 
3  and  4.  A  similar  situation  to  Fig.  6  is  seen  in  Figs.  7(a)  and 
7(b),  although  only  the  first  detector  ring  was  sensitive  to  a 
small  (<  1  mm  thick)  fat  layer  on  the  surface,  while  the  outer 
rings  sampled  deeper  and  more  radially  distal  tissue  volumes.  It 
should  be  noted  that  in  the  fixation  of  samples,  the  fat  regions 
tend  to  shrink,  "6  so  the  measurement  surface  of  these  specimens 
were  likely  flatter  during  signal  acquisition.  Also,  all  specimens 
were  cut  after  fixation  and  before  sectioning  to  make  a  given 
section  contain  only  the  interrogated  tissue  region,  so  the  fibers 
were  never  placed  over  the  edge  of  any  sample. 

The  opposite  situation  of  the  above  two  samples  is  shown 
in  Figs.  7(c)  and  7(d),  where  the  source  fiber  was  placed  over  a 
tumor  region  ~  1  to  1 .5  mm  thick,  with  normal  tissue  underneath; 
outer  detector  rings  were  placed  over  a  considerably  thinner 
tumor  layer  with  more  underlying  normal  tissue.  Taken  with 
the  above  results,  these  panels  demonstrate  that  the  spectral 
signatures  collected  with  the  SORS  probe  vary  appropriately  as 
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a  function  of  S-D  offset  according  to  tissue  type  and  location, 
not  via  any  systematic  response.  Since  some  normal  Raman 
signatures  are  present  even  in  ring  1  spectra,  Figs.  7(c)  and  7(d) 
also  shows  that  the  presence  of  normal  tissue  under  tumor  tissue 
would  be  easier  to  detect  than  the  current  problem,  since  fat  is  a 
stronger  Raman  scatterer  than  tumor  tissue.25 

Given  these  findings  regarding  sampling  depths  and  volumes, 
the  composite  spectra  were  used  for  margin  analysis  on  intact 
breast  specimens  in  the  laboratory.  Since  the  SNR  is  approxi¬ 
mately  equal  in  all  4  rings  (see  Fig.  4),  averaging  them  provides 
information  about  the  entire  sampling  volume  in  a  single  spec¬ 
trum.  This  method  also  simplifies  the  analysis  procedure;  if 
spectra  from  individual  rings  were  used,  it  would  be  difficult  to 
determine  how  to  correlate  certain  ones  with  pathology  findings. 
For  example,  although  all  spectra  in  Figs.  6(b)  and  7(b)  were 
from  tissue  sites  that  would  be  deemed  positive  margins  within 
the  spatial  extent  of  the  probe,  it  is  unlikely  that  the  innermost 
rings  were  actually  detecting  any  signal  from  tumor  tissues.  A 
possible  approach  for  using  the  individual  spectra  would  be  to  la¬ 
bel  a  measurement  site  positive  if  any  spectrum  from  the  4  rings 
is  predicted  to  be  from  a  positive  margin,  but  the  aforementioned 
correlation  issue  arises  in  the  training  of  such  an  algorithm  for 
a  retrospective  analysis.  Many  normal  looking  spectra,  like  ring 
1  from  Fig.  6(b),  would  be  labeled  as  a  tumor  and  would  likely 
cause  difficulties  for  discrimination  algorithms  trying  to  create 
decision  boundaries  between  negative  and  positive  margins. 

A  binary  diagnostic  algorithm  simpler  than  SMLR  may  seem 
like  a  more  appropriate  approach  in  this  analysis,  but  the  SMLR 
algorithm  was  able  to  significantly  reduce  the  dimensionality  of 
the  data  from  the  initial  size  of  232  variables  (one  per  3.5  cm-1 
bin)  to  perform  its  classification.  In  addition,  SMLR  provides  a 
probability  of  class  membership  that  would  be  very  useful  in  a 
clinical  application.  A  surgeon  could  act  differently  if  the  prob¬ 
ability  of  a  margin  being  negative  is  99%  versus  51%,  although 
in  either  case,  the  diagnosis  would  be  negative. 

The  results  from  using  SMLR  to  classify  the  composite 
SORS  spectra  according  to  margin  status  are  shown  in  Table 

1.  With  only  one  false  negative,  the  sensitivity,  specificity,  NPV, 
and  PPV  were  all  at  least  94%.  For  this  clinical  application, 
perhaps  the  most  important  variable  for  long  term  studies  is 
NPV,  since  a  surgeon  needs  to  be  confident  in  any  diagnosis 
of  negative  margin  status  to  prevent  recurrence  of  the  disease 
or  unnecessary  second  operations.  For  the  single  false  negative 
result  in  this  study,  the  normal  layer  overlying  the  tumor  was 
found  to  be  ~1.5  mm  thick  upon  histological  examination,  but 
prior  to  formalin  fixation,  this  layer  was  likely  around  or  slightly 
greater  than  2-mm  thick,36  which  would  surpass  the  sampling 
capabilities  of  the  SORS  probe.  It  may  also  be  possible  that  a 
slight  misalignment  between  the  probe  and  the  point  of  the  histo¬ 
logical  section  led  to  an  error  in  the  margin  size  determination. 
In  addition,  there  is  not  a  universal  standard  among  hospitals 
of  a  minimum  margin  size  required  during  breast  conserving 
surgery;  rather,  some  locations  use  2  mm,  some  use  1  mm,  and 
others  simply  require  that  no  cancer  cells  be  found  on  the  sur¬ 
face  of  the  specimen.3  We  used  2  mm  as  the  cutoff  in  this  study 
because  that  value  provides  the  best  prognosis  for  patients3  and 
is  the  most  stringent  standard  for  proving  the  value  of  SORS. 

The  classification  results  above  compare  extremely  favorably 
with  current  intraoperative  margin  evaluation  techniques.5-8  For 
example,  the  reported  sensitivity  of  “touch  prep”  is  as  low  as 


8%  (Ref.  5);  a  simple  visual  examination  has  sensitivity  and 
specificity  of  approximately  50%  and  72%,  respectively;7  and 
frozen  section  pathology,  though  its  sensitivity  and  specificity 
per  slide  are  generally  >  90%,  suffers  from  a  sampling  error  that 
brings  per-specimen  classification  accuracy  (i.e.,  overall  desig¬ 
nation  of  whether  a  second  operation  is  required)  below  85% 
(Ref.  8).  Another  optical  approach  for  intraoperative  margin 
evaluation  is  to  image  an  entire  margin  (i.e.,  one  of  the  six  facets 
of  the  “cuboidal”  excised  specimen)  at  once  with  autofluores¬ 
cence  and/or  diffuse  reflectance  modalities.  While  our  group 
has  published  a  small  study  on  the  topic,13  it  has  been  exten¬ 
sively  researched  in  recent  years  by  the  Ramanujam  group.9-12 
Using  extracted  optical  properties  of  the  tissue  from  visible 
diffuse  reflectance,  that  group  has  achieved  79%  sensitivity 
and  67%  specificity  for  discriminating  normal  from  positive  or 
close  (<2  mm)  margins  for  a  set  of  48  patients.11  Bet¬ 
ter  classification — 100%  sensitivity  and  82%  specificity — was 
achieved  by  Nguyen  et  al.  using  OCT,14  though  the  sample  size 
of  20  patients  was  much  smaller,  and  the  technology  in  its  cur¬ 
rent  state  still  relies  on  a  subjective  analysis  of  the  images.  The 
biggest  challenge  for  the  SORS  approach  presented  here,  espe¬ 
cially  compared  with  the  above  two  techniques,  is  adapting  the 
probe  and  other  system  components  to  interrogate  larger  areas 
of  tissue  in  a  shorter  time. 

The  various  optical  approaches  to  an  intraoperative  margin 
evaluation  all  hold  significant  potential  for  improving  the  stan¬ 
dard  of  care,  though  each  currently  has  its  own  strengths  and 
weaknesses.  To  date,  no  method  has  demonstrated  the  com¬ 
bination  of  sampling  speed,  volume,  and  diagnostic  accuracy 
needed  for  widespread  clinical  implementation.  The  initial  work 
presented  here  has  demonstrated  the  feasibility  and  promise  of 
using  SORS  to  evaluate  the  margin  status  on  intact  breast  speci¬ 
mens  in  a  laboratory  setting.  Studies  are  currently  underway  on 
using  the  same  approach  in  a  clinical  setting,  and  initial  results 
are  equally  as  promising  as  the  laboratory  measurements.  These 
clinical  SORS  measurements  for  breast  tumor  surgical  margin 
evaluation  will  be  the  subject  of  future  manuscripts. 
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